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We investigate theoretically the properties of an ideal 
trapped gas in a time-dependent harmonic potential. Using 
a scaling formalism, we are able to present simple analytical 
results for two important classes of experiments: free expan- 
sion of the gas upon release of the trap; and the response of 
the gas to a harmonic modulation of the trapping potential 
is investigated. We present specific results relevant to current 
experiments on trapped Fermions. 



I. INTRODUCTION 

Recently, impressive experimental results concerning 
the quantum degenerate regime of a dilute gas of trapped 
Fermionic atoms have been presented |^]. By cooling 
~ 7 X 10^ '^^K atoms to 0.5 of the Fermi tempera- 
ture, the onset of Fermi degeneracy was observed in 
the thermodynamic and in the scattering properties of 
the gas. Also, progress towards achieving fermi de- 
generacy for trapped ®Li j2| has been presented. Such 
a weakly-interacting, degenerate Fermi-Dirac gas pro- 
vides a new platform for exploring fundamental quantum 
many-body physics. Several theoretical results dealing 
with the equilibrium properties of such a gas have been 
presented For example, it has been shown that 

a two-component gas of spin-polarized ^Li undergoes a 
Bardeen-Cooper-Schrieffer (BCS) transition to superflu- 
idity at experimentally obtainable densities and temper- 
atures 

One class of experiments that is compatible with trap- 
ping protocols, and which has gleaned valuable informa- 
tion on the dynamics of dilute Bose-Einstcin condensates 
(BECs), involves monitoring the response of the gas to 
a change in trapping potential. For example, transient 
modulation of the trapping potential induces free ring- 
ing of the gas, which in the case of BECs led to a direct 
determination of low-lying regions of the quasi-particle 
spectrum . Complete release of the trapping potential 
enables one to view the free expansion of the gas; this 
established early on the essential validity of the time- 
dependent Gross-Pitaevski equation for describing the 
dynamics of latest generation of BECs which was 
later put to stringent quantitative tests [ |lO| . 

A trapped, single-component 

gas of ultracold Fermionic atoms can, for experiments 
of current interest, be considered to be ideal (noninter- 
acting) as a reasonable approximation, since atomic col- 
lisions are strongly suppressed pl| , p^ . In this paper, we 
examine theoretically the dynamics of such an ideal gas 
in a time-dependent trap. Using a scaling formalism simi- 



lar to that which has been successfully applied to trapped 
BECs [p"3[jT^] , we derive analytical results for the class of 
experiments mentioned above. The essential Fermionic 
character of the systems with which we are concerned 
is established by the initial distribution of particles in 
the trap; the time evolution of this distribution, under 
changes of the trapping potential that preserve its har- 
monicity, are rigorously equivalent to that of an ensemble 
of noninteracting particles, independent of statistics. We 
present a simple formula that describes the free expansion 
of such an ideal gas, and our results suggest a new ap- 
proach to the problem of quantitative thermometry in the 
nanokelvin regime. Due to their weak pair interactions, 
single-component Fermi-Dirac gases are attractive can- 
didates for nanokelvin thermometry; with an improved 
theoretical understanding of finite-temperature proper- 
ties of BECs , which are much more robust candidates 
for experiment at present, we can envisage a direct com- 
parison of temperatures of ultracold Bose-Einstein and 
Fermi-Dirac gases. We also examine in detail the (non- 
linear) response of the gas to a harmonic oscillation of 
the trapping frequency. This reveals a domain of driving 
frequencies and amplitudes which generates a resonant 
response of the gas. 

The non-interacting limit considered here should 
provide a good approximation for degenerate, single- 
component cold Fermi-Dirac gases. Since interactions 
in general play a smaller role for Fermi-Dirac vs. Bose- 
Einstein systems, the present approach may also provide 
a useful starting point for consideration of the dynamics 
of multiple-component Fermi gases. 



II. FORMALISM 

We start with a derivation of the equations describ- 
ing the scaling properties of an ideal gas trapped in a 
time-dependent harmonic potential. Consider a classi- 
cal particle of mass m trapped in a potential V{r,t) — 
m'^j LL!j{t)'^r'j /2 with j = x,y,z denoting the 3 spatial 
dimensions. We set ojj{t) = ojQj for t < 0, i.e. the 
trap potential is constant prior to i = 0. Newton's law 
is then expressed by pj = —mujj{t)'^rj and mij = kj. 
Invoking a scaling transformation qj = rj/"fj(t) and 
Pj = jj{t)kj — ^j(t)mrj, we obtain dr^Pj — —mw^^q^, 

with Tj{t) = J* dt''-fj{t')^^ , if the scaling parameters 7j(t) 
satisfy the equations 
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In quantum mechanics, the Heisenberg equation for a 
non-interacting gas in a time-dependent harmonic po- 
tential takes the form: 



2m 



W^ + V{r,t)]i;{r,t). 



where ip{r) is the field operator for an atom in a sin- 
gle hyperfine state at position r, which obeys the usual 
Fermi-Dirac anticommutation relations. The quantum 
mechanical analogue of the rescaling outlined for the clas- 
sical case above corresponds to writing |13 



(3) 



with qj{t) — rj/jj{t), as defined above. If each -fj satis- 
fies Eq.(|l|), by writing $[q(i)] = Tlj4'[qj{t)], we obtain 



2m 



-mLjl^q^]^[q^{t)]. (4) 



Thus, the time-dependent problem has been reduced to 
the trivial case of evolution in a time-independent har- 
monic trap in the rescaled variables (t^ ,qj). Determi- 
nation of these variables requires only the solution of 
the three ordinary differential equations expressed by 
Eq.(|l]); the three spatial dimensions can be treated in- 
dependently. An appealing qualitative description of the 
solution of these equations, cast in the context of one- 
dimensional scattering theory, has been given by Kagan 
et al. From Eq.(||), it follows that the density p(r, t) 
of the gas for a given time t is given by 



p(r,i)= (V'Ur,t)V'(r,t)) = 



1 



Po[q(t)] (5) 



where po{r) is the particle density for t = 0. Thus, we 
can calculate p(r, t) for any time t > for modulated fre- 
quencies iOj(t), by solving Eq.(^, subject to the bound- 
ary conditions 7j(0) = 1 and 7^ — 0. Equation (|l]) also 
describes the two-dimensional, non-ideal BEC subject to 
isotropic variations of the trap potential, with Eq.(@) also 
being applicable if po is obtained by a solution of the 
Gross-Pitaevski equation jrif . 

We now analyze the solution of Eq.(|l|) for two cases 
of experimental relevance: free expansion; and harmonic 
modulation of the trapping potential. 



III. FREE EXPANSION 

To model a free expansion experiment, we take 0Jj{t) — 
for t > 0. Solving Eq.(0) with the requirement that the 
gas is in equilibrium for f < 0, we obtain 7j(i) = 1 for 
t < and 



7,(t) 



1 



t 



(6) 



for t > 0. The time-dependent width of the cloud af- 
ter the trap has been dropped is given by ^ {rj){t) = 

7j(t)w (f|)(0). To describe the aspect ratio a{t) of a 



(2) cylindrically-symmetric cloud, we find 



a{t) 



(z2)(t) 



«(0)i 



t2 



1 + ujLz' 



(7) 



We now apply these results to an ideal gas in two limiting 
cases. 

We first treat the case of T = 0, with the chemical 
potential fJ^piT — 0)/hujQj ^ 1 for j = x,y,z. This cor- 
responds to the semiclassical limit, appropriate to cur- 
rent experiments with more than a few hundred trapped 
atoms. In this limit, the initial density profile is well 
described by the Thomas- Fermi (TF) approximation 
This gives the integrated density p[x,z,t) = J dyp{v,t), 



p{x,z,t) 



Anh UJoylxJz 



(8) 



with = woz/woa: and Rp ~ \j2pp jrauJ^.^. We see 
that the cloud becomes isotropic for woji ^ 1- This 
is the result one would get for a classical gas at fi- 
nite temperature, in which the momentum distribution is 
isotropic, and is consistent with the initial isotropic mo- 
mentum distribution implicit in the TF approximation 
(in BECs in anisotropic traps, on the other hand, the 
long-time expansion is anisotropic due to macroscopic 
population of an anisotropic initial state, a key result in 
their first observation ||l8|). The initial aspect ratio of 
a(0) — \z and Eq.(|^) evolves to Q!(i) ^ 1 as i ^ 00. In 
Fig.(|l|), we plot the integrated density p{x, z, t) for T = 0, 
Hp{T =^0) ^ 20;iwoa;, and t = (a) and t ^ 20/luox (b). 
We have taken A2 = 19.5/137 and Aj, = ujoy/ujox — 1 
corresponding to current experiments on trapped '^''K 
atoms 
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FIG. 1. Contour images of the integrated density p{x, z, t) 
for a freely expanding gas. There are 9310 atoms in the cloud. 
We have defined = {h/muiox)^^^ . 

The two snapshots of the integrated density clearly 
show that the gas becomes isotropic for large expansion 
times. This is further illustrated in Fig.(^), where we plot 
a(t) as given in Eq.(^ for the same set of parameters. 
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FIG. 2. The aspect ratio for a freely expanding gas. The 
solid line is for = 19.5/137, = 1 and \ip — 20fttJoa;- The 
dashed line is for Az = Aj, = 1/50 and fiF = huox. 

Note that a{0) = Xz for any initial density distribution 
of the form po{r) = fiJ2j^oj''^j)- Such initial distribu- 
tions will become isotropic in an free expansion experi- 
ment. For instance, for T > Tp = ^f{T = Q)/kB, the 
density is well described by a classical gaussian profile, 
i.e. po(r) « exp[-P{fiF-mY,. Lul^r"^ /2)] andFig.(|) thus 
describes a trapped gas of fermions for any T within the 
TF approximation. This means that one cannot detect 
the onset of Fermi degeneracy by looking at the aspect 
ratio of the expanding gas. 

As an example where po{r) 7^ '^oj'^i )' ^"^^ 

consider a case with pp/huJox < 3/2 such that only one 
level is occupied in the x-direction and ^ip/^^oj ^ 1 foi' 
j — y, z. The gas is initially strongly confined in the 
x-direction. The integrated density profile is then 

Pi..y,i)o.e^^nl.llr^pflAY (9) 

\^ 2^_F - Tiujox j 

for t > 0, yielding a{0) — \/ZhiOQxXz j \/2[if — fiuJox- In 
Fig.(^), we plot the aspect ratio for a free expansion for 
p,F — hujox and Aj, = Az = 1/50 using Eq.(^. We have 
a(0) = \/3/50 and a{t) ~~> \/3 for t co. The gas, 
which initially is strongly confined in the x-direction will 
for ujajt 3> 1 become most confined in the y-, and z- 
directions. This is, of course, a direct reflection of the 
uncertainty relation giving higher average momentum in 
the a:-direction. However, observation of such a quantum 
effect requires a highly anisotropic trap. For the case of 
lO** trapped atoms we would require Aj, = A^ ~ 1/250. 
Thus, for Fermi-Dirac particles the anisotropy of the ex- 
panded cloud due to the Heisenberg uncertainty principle 
is considerably less than that encountered in the Bose- 
Einstein case. 

From a measurement of the density at any time t un- 
der free expansion it is straightforward to determine the 
initial density distribution pq{v) using Eq.(||)-(|6|). This 
suggests that a trapped, single-component Fermi-Dirac 
gas could serve as a low T thermometer. Assuming the 
gas is in thermodynamic equilibrium for t < 0, one could 
infer the temperature of the gas from a knowledge of its 
density distribution, the calculation of which is a trivial 
problem of summing over fractionally occupied trap lev- 
els. Hence, a determination of T is simply a matter of 
fitting a calculated density to the measured p(r) . Density 



profiles for an ideal Fermi-Dirac gas for various tempera- 
tures have already been given elsewhere However, if 
only one hyperfine level is present at all times, the atomic 
collisions required for thermalization are suppressed at 
low temperatures, and evaporative cooling of the gas will 
be inhibited. This can be remedied by initially trapping 
the gas in two hyperfine states, to allow s-wave collisions 
between atoms in different states. One can then deplete 
the trap of one hyperfine state by applying a microwave 
field that depletes one of the states. If this depletion 
occurs on a time scale that is long compared to colli- 
sion times in the sample, the remaining atoms will be in 
thermodynamic equilibrium, and their expansion can be 
described by the results in this section. In fact, a proce- 
dure of this type is being used in present experiments on 
trapped atoms ||l|. 

Alternatively, one could drop the trap with both hyper- 
fine levels present. The initial density profile and the free 
expansion of the gas will then be perturbed away from the 
ideal gas result by the interactions. When the effect of in- 
teractions is small, the initial equilibrium density and the 
subsequent expansion can be calculated relatively easily. 
A calculation of the effect of the interactions on the equi- 
librium distribution has been presented elsewhere, Q and 
the effect of the interactions on the free expansion can be 
calculated within the mean-field approximation, starting 
from the simple uncoupled solutions given by Eq.(^ for 
each spatial direction |19|] . 

IV. DRIVEN OSCILLATIONS OF THE GAS 

As noted above, modulation of the trapping frequency 
has proven to be a useful technique for probing the low- 
lying collective modes of trapped BECs. We therefore ex- 
amine theoretically the non-linear response of a trapped 
gas to such a modulation. The non-interacting case 
treated here is in some sense opposite to the hydrody- 
namic limit we have treated elsewhere ||20| . To model a 
typical experiment, we assume that the trapping frequen- 
cies take the form LOj{tY = ij^ojll ~ 2rj cob{u) nt)] for i > 
with being the driving frequency and rj the driving 
amplitude. Instead of solving Eq.(|lJ) with this form for 
t-^jit), it turns out to be easier to go back to the origi- 
nal classical equation of motion rj — —Ljj{t)'^rj by using 
rj{t) = qj{t)jj{t) and qj{t) = exp[±iuJojT{t)]. Thus, by 
writing ^j(i) — 7_; (i) exp[ia;ojT(<)], and x = ^nt/^, we 
obtain: 

dli, + [a-2qco^(2xMj=0 (10) 

with a ~ 4/a)^, q = Arj/uP and uj — LOg/bJoj. Equation 
( |lO| ) is a variant of Mathieu's equation [gl| , whose prop- 
erties have been extensively studied. Its behavior is easy 
to explain in one case of experimental interest: that when 
the trap frequency is modulated during a finite interval 
< t < to, and then returned to its original value. The 
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subsequent motion of the cloud is described by a solu- 
tion of the time- independent problem LUj (t) = ujQj with 
arbitrary initial conditions: 



7iW = \/ VE^ ~ 1 sm{2ujojt + c)+E, (11) 

where E = {ujQ^^jj+-fj+'j^'^)/2 > 1 is a conserved quan- 
tity for i > to; c = arcsin[(7gj- —E)/y/E^ — 1]; and 7oj = 
Ijiio) is the value of ^j{t) immediately after the mod- 
ulation ceases. Equation ( pi] ) has a discrete frequency 
spectrum (frequencies of 2nuJoj with n = 0,1,2...), as 
expected for a non-interacting gas in a harmonic trap. 
The linear response limit is recovered by letting i? — s- 1+ 
inEq.(|ll|), to obtain 7j(i) = 1-1- S sm{2ujojt). 

In essence, the problem of predicting the response of 
the gas to a harmonic driving with frequency cj/^ and 
amplitude 77 is reduced to an analysis of the well-known 
solutions to Mathieu's equation. In the parameter space 
(a,(jr), there are regions where the solutions of Eq.(p^) 
are unstable, i.e. their amplitude increases exponentially 
with time. Also, there are stable regions where the solu- 
tions remain bounded. The solutions on the boundaries 
between these regions are the Mathieu functions . Us- 
ing |7j(t)| = \£_j{t)\, a = 4:/uj'^, and q = ^rj/u}"^, this means 
that for certain regions in the (wd, ?7)-space, the response 
of the gas diverges as the driving time to increases, i.e. 
there is a resonant response, whereas in other regions 
the response of the gas remains finite for any value of Id . 
Specifically, for a — n? with n = 0, 1, 2 . . ., the solutions 
to Eq.(p^) diverge in time for an arbitrary small q pl| . 
Thus, for 



ujd — "iuJoj /n, n=l,2,3. 



(12) 



the response of the gas is resonant for an arbitrarely small 
driving amplitude and the amplitude of its oscillations 
will diverge with the driving time. In terms of Eq.(ll), 



if the gas remains trapped in the time-independent po- 
tential LiJj{t) = LOQj after the driving [t > i^), we have 
E 1 and the oscillations of the gas will be large and 
contain many harmonics. The resonance for ojd = 2a;oj 
is, of course, the usual excitation frequency for an even- 
parity perturbation for a non-interacting gas. However, 
the resonances for = woj, 2ti;oj/3, woj/2 etc. do not 
correspond to new modes. They simply reflect the fact 
that for these driving frequencies, the trapping potential 
is doing resonant work on the gas. Note that the present 
scaling formalism does not predict the higher-frequency 
modes at n2u)Qi with n > 2, as they do not correspond 
to simple dilations/contractions of the cloud. 

We will now examine the width of the unstable (reso- 
nant) regions around lud — 2uj(),j/n for 77 — > 0. For the 
ljd = 2ccioj resonance, the unstable region of the Mathieu 
equation is bounded byl — g<a<l-|-(7 for g ^ . 
This gives, that for driving frequencies and amplitudes 
such that 



the response of the gas is divergent. Hence, the resonance 
region of the gas for finite but small driving amplitude 
77 has a reasonable width and should be relatively easy 
to access experimentally. Likewise, for the ljd = Lo^j 
resonance, the unstable region is bounded by 



1 - 5?7^/6 < a) < 1 ?7^/6, 77 < 1 



(14) 



and it should be experimentally accessible. For the lower 
frequency resonances (tj/j = 2ujQj /n, n > 3), it turns out 
that the resonance regions are very narrow for 77 ^ 
as the boundary lines between the stable and unstable 
solutions of Eq.(p^) only differ by terms of order or 
higher |^l|. Thus, these resonances are difficult to find 
experimentally for 77 —s- 0. Instead, one could increase the 
driving amplitude 77 for a given frequency ljd] for large 
enough 77, an unstable region is reached and the response 
of the gas diverges. The above results are illustrated 
in Fig.(||), where we plot the response of the gas as a 
function of driving frequency lod and amplitude 77. 




2- 



7/ < a) < 2 + 77, 7/<l, 



(13) 



FIG. 3. The regions of stability/unstability for the re- 
sponse of a non-interacting gas to a modulation of the trap- 
ping potential with frequency ljd and amplitude r]. 

The lines in the plot separate regions where the re- 
sponse of the gas remains finite from regions where the 
oscillations of the gas diverge with driving time. For ex- 
ample, if one modulates the external potential with an 
off- resonance frequency of a) = 5/3, the response of the 
gas remains finite for small 77. However, if one gradu- 
ally increases the driving amplitude 77 for this frequency, 
then one should observe a divergent response of the gas 
for 77 larger than ^ 0.33. The plot is generated using 
the well-known properties of the Mathieu equation and 
the mapping {a,q) (4/(1;^, 4r//(i^). Note the unstable 
(resonance) regions for r/ ^ for lojj = 2ujQj/n. The 
regions have zero width for 77 ^ for 77, > 3 as pre- 
dicted above. For increasing driving amplitude 77, the 
unstable regions grow as expected. For lod —^ 0, the 
response of the gas for rj <^ 0.5 is finite whereas it di- 
verges for f] >~ 0.5. Physically, this corresponds to the 
fact that for ujD/uJoi <^ 1 and 77 < 0.5, the gas returns 
adiabatically to its initial state after one period of trap 
modulation and there is no net work done. Contrary, for 
77 > 0.5 the trapping potential becomes inverted for cer- 
tain times resulting in a divergent response of the gas. 
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Mathematically, the transition region for 77 ~ 0.5 conies 
from the fact, that for a ~ 2^, the even and odd periodic 
solutions (Mathieu functions) of Eq.([lO|) start to differ 
in "energy" (the Mathieu characteristic value a HJ]), as 
the tunneling between successive minima of the potential 
potential cos(2x) becomes significant. 

Again, since a spin-polarized gas of cold fermionic 
atoms is effectively ideal, the results presented here 
should apply to the response of such a gas to a harmonic 
modulation of the trapping potential. 



V. CONCLUSION 

In conclusion, using a scaling formalism we have been 
able to derive simple analytical results for the dynam- 
ics of ideal gases trapped in time-dependent harmonic 
traps. We have concentrated on two important classes of 
experiments: For free expansion, we show how the initial 
density profile of the gas can easily be determined from 
a measurement of the density profile at any time t after 
the trap has been dropped. We proposed a low T ther- 
mometer based on these results. Also, we showed how the 
problem of determining non-linear response of the gas to 
a harmonic modulation of the trapping frequency can be 
mapped to an analysis of the well-known properties of the 
solutions to Mathieu's equation. We identified regions in 
77)-space where the response of the gas was diver- 
gent with the modulation time. Especially, we were able 
to predict unstable regions for 7/ — > reflecting the fact, 
that the trap is doing resonant work on the gas. Since 
ultracold spin-polarized fermions are non-interacting to a 
very good approximation, our results should be directly 
relevant for current experiments in this very active field 
of research. 
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